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Abstract 

We discuss the quantum-classical correspondence in a specific dissipative chaotic system, Duffing 
oscillator. We quantize it on the basis of quantum state diffusion (QSD) which is a certain formu- 
lation for open quantum systems and an effective tool for analyzing complex problems numerically. 
We consider a sensitivity to initial conditions, " pseudo-Lyapunov exponent ", and investigate it in 
detail, varying Planck constant effectively. We show that in a dissipative system there exists a cer- 
tain critical stage in which the crossover from classical to quantum behavior occurs. Furthermore, 
we show that an effect of dissipation suppresses the occurrence of chaos in the quantum region, 
while it, combined with the periodic external force, plays a crucial role in the chaotic behaviors of 
classical system. 
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I. INTRODUCTION 



The quantum-classical correspondence is a very important problem related to the foun- 
dation of quantum mechanics. However, it is difficult to consider this correspondence if 
a classical system which we should quantize is a chaotic system (or generic nonintegrable 
system). For, the phenomena in the quantum mechanics corresponding to the chaos in the 
classical mechanics can not exist generally due to quantum mechanical dynamics Qj. 

Nevertheless, the various works have been done in Hamiltonian systems. The central task 
for studying the quantum-classical correspondence in Hamiltonian chaotic systems is how 
the measured quantities in a quantum system relate to the information of trajectories in its 
corresponding classical system. The study in this point is very fruitful P, 0, 1^1 . 

However, researches limited in Hamiltonian systems are not sufficient to study the 
quantum-classical correspondence in chaotic systems. There exist another types of chaos 
not having Hamiltonian. This phenomenon can occur in dissipative systems that a definite 
Hamiltonian does not exist. In this paper, we discuss the quantum-classical correspondence 
for Duffing oscillator in the view of an open quantum system. 

There are several reasons why we study this system in order to consider the dissipative 
quantum chaos. The authors in Ref. Jd] reported interesting results for the same system 
as our model. Especially, the proof of existence of chaotic behavior in their classical limits 
is very important. We think, however, their results are insufficient to discuss the crossover 
from classical to quantum behavior. In fact, we show that a new method, an analysis based 
on "pseudo-Lyapunov exponent", is possible to clarify the crossover behaviors. Moreover, 
there are experimental propositions for dissipative quantum anharmonic oscillator which are 
not identical with the Duffing oscillator. In Ref. , the analysis is worked in such a model. 
Furthermore, the investigation of this problem is related to not only the quantum chaos but 
also the several fundamental problems in the quantum mechanics, for example, the influence 

Finally, it also has interesting features in the classical mechanics [lOj. The other svstem for 
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the dissipative quantum chaos, e.g., the dissipative kicked top, is studied in Ref. 

The main aim is what happens in the stage between quantum and classical regions in a 
dissipative system, as the Planck constant changes effectively. It is expected that the various 
phenomena related to quantum-classical correspondence should occur. But it is not clear 
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what kind of quantity we should investigate to find such a crossover. We examine a quantity 
sensitive to the initial condition and define "pseudo-Lyapunov exponent" as its candidate. 
Then we will discuss such a crossover in the quantized Duffing oscillator. The effective 
Planck constant j3 and the effective Planck cell play an important role in this analysis. 

We obtain several fruitful results. First, we find that "pseudo-Lyapunov exponent" is 
certainly positive when the system of Duffing oscillator is in a classical region. Furthermore, 
we show that there is a certain clear critical stage in which the crossover from classical to 
quantum behavior occurs. We find also that the effect of dissipation is different between 
classical and quantum regions; this, together with a periodic external force, plays the essen- 
tial role for chaotic behaviors in the classical region but suppresses the occurrence of chaos 
in a quantum region. 

The contents of this paper are as follows. In Sec. m we explain the model. Duffing os- 
cillator. Especially, we explain the phenomenological description of open quantum systems, 
quantum state diffusion (QSD). Our method is identical with that in Ref. j^. We explain 
in detail the introduction of scaling parameter (3 which is very important for the investiga- 
tion of crossover behavior. In Sec. IIIH we show several numerical results without averaging 
over the ensemble for complex Wiener process used in QSD. These methods give the proper 
results around /? = 0.01. Such a reserch has been already studied in Ref. p. The results 
are not first realizations, but are important to explain our motivation of the analysis in the 
next section. In Sec. IIV| we show the main results in this paper. We introduce a quantity 
sensitive to initial conditions, and find that there exists a clear crossover from classical to 
quantum behavior as /5 — > 1. In Sec. El we summarize this paper. In Appendix El we derive 
an equation used in Sec. IIVDI 



II. MODEL 



In this section, we first review the classical Duffing oscillator briefiy. We quantize the 
Duffing oscillator as an open quantum system phenomonologically, using QSD. This method 
is identical with that in Ref. jy]. 

The equation of motion for classical Duffing oscillator is the following: 

X + 2Tx + - X = gcos{nt). (1) 



3 



It is known that the various behaviors can occur depending on the set of parameters F, g 
and VL. The chaotic behavior appears in the case of r = 0.125, g = 0.3 and Q = l.OoQ- We 
find a strange attractor in Poincare surface, which is obtained by putting (x, p) in a phase 
space by every interval of 27i/Q. Hereafter we use this set of parameters. The appearance 
of strange attractor in such a surface is one of the properties in dissipative chaotic systems. 

We describe phenomenologically the dynamics of such a system without a well-defined 
Hamiltonian; we regard it as an open quantum system. We assume the Markovian dynamics 



and choose the effect of dissi pat ion phenomenologically 
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Then the dynamics of system 



is described by QSD jy, [IJ, [l^ [ly, [TO, llSl in which the pure state vector of system 
evolves according to the Ito stochastic differential equation: 

ih y 2 

-^{L^){L)ym+[L-{L)y^)dt (2) 

where M{d^} = 0, M{d^d^} = and M{d$,*dC,} = dt. d^ describes the increment of complex 
Wiener process and M expresses the ensemble average for it. (•) represents {ip\ • \ip). The 
quantum expectation value for an operator O is represented by 

M{(V'|0|V')} = Tr{dp}, 

where p is a reduced density matrix for the system 3^. ^ is a certain self-adjoint operator 
and it is called by Hamiltonian since the term related to H describes the unitary evolution. 
L is called by a Lindblad operator and describes the effect of dissipation. The QSD is 



equivalent to the Lindblad master equation [20, |2J| : 

p = -'^[H,p] + LpD - h^Lp - ^pDL (3) 

The Lindblad master equation is a quite general formulation for open quatum systems 
satisfying the Markovian dynamics, trace preserving and complete positivity. The QSD is 
a very effective tool for numerical simulation of complex problems [6], compared with the 
description depending on the master equation. In this paper, we use the algorithm for QSD 
which the authors in Ref. invented. The QSD is also possible to explain a measurement 
processes in the quantum mechanics 3, • 

In order to describe the dynamics of Duffing oscillator in the quantum mechanics, we 
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define the Hamiltonian H and a Lindblad operator L in tlie Eq. (j2I) as the followings: 

H = Hd + Hr + H,,, (4a) 

Hb, = -{xp + px), (4c) 

Hex = —gmlujQX cos^ut), (4d) 



x + -p. (4ej 



The Eq. (Pc|) means the strength renormalization for coupling of interaction between system 
and environment. In fact, if we implement the canonical transformation, x ^ x and p — >■ 
p — 7mx, we can obtain the following equations: 

^^2 , _ (. , T^A ^2 



\ cjo / V muJoh 
The Eq. fl4d|) means the external force depending on time periodically. Notice that the right 

hand side of Eq. Q is independent of the time. Therefore, it is a difficult problem how the 

generator depends on the time. We determine simply it so as reproduce the external force 

in the equation of expectation values. 

We rewrite Eq. (j2|) into the dimensionless form: 

-^{k^){k)yij)dT+ (^k - {k}y^|J)dc (5) 

where M{dC} = M{dCdC} = 0, M{dCdC} = dr, Hp = H/fkuo, k = L/^, r = uj^t 
and dC, = ^/uJod^. We define the unit of energy as hjjQ. Moreover, we redefine x and p as 



Q = a/ mujQ I h X and P = ^/TJmw^ p, respectively. Thus we obtain the dimensionless 
Hamltonian Hp and Lindblad operator K: 

HfS = HD + HR + Hex, (6a) 
Hd = \P' + ^Q' - \q\ (6b) 

Hr=^[qP + PQ), (6c) 

Hex = ~Qcos{nt), (6d) 

k = ^{Q + iP^, (6e) 
where f2 = uj/ooq, T = 'y/ujQ. The is the ratio of h to the characteristic action of system, 

We can effectively change /i, varying (3'^ in the numerical computation. Notice that the 
varying of (3 is just the scale transformation for system: AQ/MKQ)} and AP/M{(P)} 
should vanish when (3 goes to zero, where AQ^ = {Q - {.Q)f, AP^ = {P - {P)f, = 
M{(Ag2)} and AP = Jm{{AP^)} Q- 
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We consider (3=1 for a moment. Using Ito calculus, we obtain the following equations 



d{Q) = {P)dT + Vf[[(vQ-^-^ 

+iVQp}dC + c.c], (8a) 



d{P) = {-2T{P) - {F) + g cos {nT))dT 

+Vf[f^VQP + t(Vp-^^yC + C.C.] , (8b) 

where Vq = {AQ^), Vp = (AP^), 2Vqp = {AQAP + APAQ) and F = Q'^ -Q. Notice 
that, if we can approximately neglect moments more than second order, then we can see 
that Eq. (jSI) reproduces the equation of motion for classical Duffing oscillator. This is one 
of physically useful advances in the QSD. In such specific realization of stochastic 

process does not allow any deviation from the classical behavior. Therefore, we are able 
to guess reasonably that the classical regions would be robust for a specific realization of 
stochastic process. 



III. BEHAVIOR IN THE CLASSICAL REGION 



The discussion in Sec. HT] allows us to obtain the proper results around P = 0.01 without 
averaging over the ensemble for complex Wiener process C (''")• See, Eq. (jSj). In this section, 
we show the several numerical results without averaging over C(^)- These results are useful 
to understand the behavior of quantized Duffing oscillator intuitively for different values of 
(3. 
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FIG. 1: These are the stroboscopic maps for {{Q), {P))- The each point in these figures represents 
the data at every 2Tr/Q for a single realization of complex Wiener process. Figures (a), (b), (c) and 
(d) are for /? = 0.01, 0.10, 0.40 and 1.00, respectively 



First, we show the stroboscopic maps for {{Q), (P)) in the Fig.^for a certain realization of 
C{t). The each point in these figures represents the data at every 27c/Q for a single realization 
of ({t). The initial state is a pure coherent state \a = 0){a = 0|, where Re{a} = V^{Q) 
and Im{a} = ^/2{P). These show that a strange attractor appears certainly and the system 
behaves chaotically in /? = 0.01, while it has been lost in /? ~ For intermediate case, 

there remains the remnant of strange attractor. We find that the scale of system gets large 
as P goes to zero. These observations are successful to show the loss of chaotic behavior 
except for j3 = 0.01 at least. Therefore, let us call that the system is in the classical region 
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for (3 = 0.01 and in the quantum region for (3 = 1.00, respectively. 
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FIG. 2: These are the stroboscopic maps for ((AQ^), (AP^)). The each point in these figures 
represents the data at every 27r/il for a single realization of complex Wiener process. Figures (a), 
(6), (c) and id) are for /? = 0.01, 0.10, 0.40 and 1.00, respectively 

Next, we show the stroboscopic maps for ((AQ^), (AP^)) in Fig. |2l with the same initial 
condition, \a = 0){a = 0|. Since we only investigate the time evolution of stochastic pure 
state in the numerical computation, we compute these quantities without averaging over 
ensemble for C(^)- These results indicate that the state almost preserve minimal uncertainty 
relation for every j3. These results, combined with Fig. also verify the argument in Sec. HTl 
numericaUy: AQ/ (Q) and AP/ (P) goes to zero if (3 goes to zero. 

These analyses without averaging over the ensemble for ({t) have been already studied 
in Ref. jsl. Especially, Fig. ^ agrees with the results in it. The results in this section are 
not first investigation using QSD. However, these are important to explain our motivation 
of the analysis based on " pseudo-Lyapunov exponent " in the next section. 
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IV. CROSSOVER FROM CLASSICAL TO QUANTUM BEHAVIORS 



In Sec. mil we have shown numerically that the quantized Duffing oscillator preserves the 
characteristic property of the classical dynamics, chaos, in the case of (3 = 0.01. However, 
this result is inadequate to understand fully the quantum-classical correspondence in this 
model due to the following points. First, it is doubt whether in j3 = 0.01 the chaotic dynamics 
survives or not, since we only obtain a figure like strange attractor in the stroboscopic maps. 
The problem remains even if one claims on the basis of this assertion that the chaotic 
behavior may occurs in /3 = 0.01. It is not clear at what region of intermediate f3 = 0.01 
and f3 = 1.00 the classical behavior survives. The definition of the classical region or the 
quantum region is obscure. Finally, we do not consider the proper quantity related to the 
behavior of system, as /3 — > 1.00. For example, {Q) and (P) in Fig. ^ are not the ensemble 
average Tt{Qp} and Tr{Pp} respectively, are not the values measured in an experiment. 

In this section, we introduce a quantity sensitive to initial conditions, "pseudo-Lyapunov 
exponent" . We examine the above three points, based on the analysis of this quantity. Using 
this, we can find that there exists the clear crossover from classical to quantum behavior as 



A. "Pseudo-Lyapunov exponent" 

We consider the quantity corresponding to instability of classical trajectories which is the 
most conspicuous characteristic in the classical chaotic systems. We define the separation 
of trajectories A(r) as the following equation: 




(9a) 




(9b) 




(9c) 




(9d) 
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where i = 1,2. Pi(t) and P2(t) denote two density matrices for different initial states 
pi(0) and ^2(0), respectively. Actually, these initial states are the pure coherent states 
|Q;j)(Q;j|, (i = 1, 2). ctj is related to initial condition (Qj, Pj) by at = \/2{Q^ + iPi). Eq. (jHal 
represents the distance in M{{Q)}-M{{P)} plain. The subscript of {1, 2} in Eq. (jHal rep- 
resents the summation over the sets of chosen initial conditions and N is the number of 
those sets. The behavior of A(t) is the sensitivity to initial conditions. We investigate this 
behavior in detail, varying j3. 

B. Effective Planck cell 

We have to choose a suitable value as the separation e = A(r = 0) of two different initial 
conditions, before the numerical simulations. Notice that two points in the phase space are 
not distinguishable in the view of quantum mechanics, if they coexist inside the same Planck 
cell. The Planck cell is limited by the Heisenberg's uncertainty relation. 

We explain what is called an the effective Planck cell. In this model, the commutator 
[Qj-P] = [x,p]/j3'^S = 1 is fulfilled. Then the Planck cell has a constant volume of AQAP = 
1/2 in the scaled phase space, whereas it has AxAp = h/2 = j3'^S/2 in the original phase 
space. With the fixed value of typical action S for the system, the smaller (3'^ corresponds 
to the smaller h and the system exhibits the more classical behavior. Thus we can define an 
effective Planck cell as (3'^S/2 ; its linear size is almost equivalent to (3 in the unit of a/ S/2. 

The concept of effective Planck cell gives us how e should be determined. We can consider 
two kinds of determination of e for (3: (1) e = 0.01 (fixed), where two points in the phase 
space are distinguishable only for the classical region (/3 = 0.01). (2) e ^ f3, where those are 
distinguishable for each (3. 

C. Crossover behavior 

In the first, we show that the results of simulation of A(r) with e fixed as 0.01. The initial 
two points are separated at a distance of effective Planck cell for (3 = 0.01 and coexist inside 
the cell for the other cases {(3 > 0.1). In Fig. Efa), we find an exponential increase of A(r), 
a characteristic behavior of chaos. This corresponds that maximal Lyapunov exponent is 
positive in classical mechanics. This behavior is also consistent to Fig. C^a). These two 
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FIG. 3: These figures are the time evolution of A(r) with e fixed as 0.01. Figures (a) and (b) are 
obtained with single realization of complex Wiener process for each initial condition (20 samples). 
Figures (a), (b), (c) and (d) are for (3 = 0.01, 0.10, 0.40 and 1.00, respectively. 

facts verify that the quantized Duffing oscillator keeps still a chaotic behavior for (3 = 0.01. 
In Fig. 121 we see very different behaviors between (b) and (c)-(d). For these values of 
/?, all initial two points coexist inside of the effective Planck cell and are indistinguishable 
from each other. Nevertheless, starting from the inside of the effective Planck cell, A(r) 
for (3 = 0.10 increases gradually and crosses the size of effective Planck cell after some 
duration and then increases simply. This suggests that the remnant of chaotic dynamics 
still survives for f3 = 0.10. However, A(r)s for f3 = 0.40 and 1.00 always stay within the 
effective Planck cell. The chaotic dynamics has been completely lost in these cases. This 
observation suggests that there exists some critical stage as P goes from 0.10 to 1.00. We 
argue that there is the crossover from classical to quantum behavior around (3 ~ 0.40 due 
to Figs. H (b)-(d) together with Figs. □ (b)-(d). 

Let us show other results in the case of e ~ /?, where the initial two points are separated 
at a distance of the effective Planck cell size. We compute A(r) for (3 = 0.10, 0.40, 0.60, 
1.00, 1.50 and 2.00, respectively. In Fig. 0] (a), it is shown that the behavior of A(r) for 
f3 = 0.10 increases exponentially, which is similar to one in Fig. |21(b). In Figs. E] (b)-(f), it 
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is shown that each A(r) takes larger values than the size of effective Planck cell only for a 
very short period. After this period, each A(r) always stays within the cell. Thus we insist 
again that there is the crossover behavior at /3 ~ 0.40. 
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FIG. 4: The time evolution of A(r) with e = /5. Aasymp is the asymptotic value of A(r). D is 
the right hand side of Eq. (|1U|) . Figure (a) is obtained with single realization of complex Wiener 
process for each initial condition (20 samples). Figures (b)-(f) are obtained with averaging over 
100 realizations of complex Wiener process for each initial condition (10 samples). Figures (a), 
(b), (c), (d), (e) and (f) are for /3 = 0.10, 0.40, 0.60, 1.00, 1.50 and 2.00, respectively 
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The asymptotic value of A. 



D. Analyses in quantum and deep quantum regions 

We investigate the case of e ~ /3 in further detaiL Especially, we focus on the discussion 
about the case of /3 > 0.40. Hereafter, we call the case of (3 > 1.00 the deep quantum region. 

In Figs. 0] (b)-(f), it is found that, except for a very short period after starting time, 
A(r) for each f] decreases for some duration and tends to approach a certain constant value 
^asymp asymptotically, which is indicated by the dashed lines. In Fig. these asymptotic 
values vs. /? are shown. Let us put down Tasymp and tq, respectively, when A(r) approaches 
^asymp and crosses down the boundary of effective Planck cell. 

Then, using the QSD, we can obtain an approximate upper bounds of A(r) for t > tq as 
following: 



In Appendix El "we derive this estimation in detail. Curves calculated by the right hand 
side of Eq. pO|l . -D(t), are represented by broken-dotted lines in Fig. |3J which show that 
Eq. (fTn|l is a good approximation of the upper bound for tq < r < Tasymp- In the derivation of 
Eq. fllO|) . two assumptions are used. First, we assume that the value of A(r) is smaller than 
the size of effective Planck cell. Secondly, we assume that the dissipative effect dominate over 
the systematic time evolution by the Hamiltonian for some time duration from the starting 
time. Thus, A(r) for tq < r < Tasymp is probably described mainly by the dissipative effect. 
Eq. (fTn|l should not work well for r > Tasymp, where if = is not a good approximation and 
there exists many other sources confining A(r) within a small value. A(r) may become the 
constant value, /S.asympi related to the inherent property of system. In classical mechanics, 
the dissipative chaotic systems like Duffing oscillator generate the chaotic dynamics due to 
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FIG. 6: The asymptotic time Tasymp- 

the coexistence of dissipative effect and periodic external force. In the quantized Duffing 
oscillator, if the action of system is much greater than h, it seems similar to classical, i.e., 
the existence of dissipation is very important for occurrence of chaotic dynamics. On the 
other side, if the action of system is smaller than h, that is, in quantum and deep quantum 
case, the analysis of Aasymp and Tasymp suggests that the effect of dissipation suppresses even 
the occurrence of chaotic behavior. 

In Fig. El it is found that Aasymp takes a certain constant value not depending on (3 
in the deep quantum region. According to the above discussion, we find that A(r) well 
characterizes the behavior of system between the classical and quantum regions, but not in 
the deep quantum region. It will be necessary to investigate other quantities, for example, 
the higher moments and Wigner function, in order to analyze the behavior in the deep 
quantum region in detail. However, A(r) is an effective quantity enough to investigate the 
crossover behavior between the classical and quantum regions. 

Finally, we explain our method estimating Tasymp- First, we determine Aasymp and tq for 
each f], using the results of simulation. Secondly, we estimate Tasymp based on the relation 

D{Tasymp - Tq) = Aasymp foT Cach (3. The rCSuhs are To, Tasymp) = (0.40, 5.71, 8.63), 

(0.60, 2.88, 11.9), (1.00, 1.45, 19.4), (1.50, 1.04, 23.5) and (2.00, 0.91, 25.9). It is found that 
Tasymp takcs a coustaut value not dependent on P in the deep quantum region. This result 
is consistent with an approximation of D at large r: 



Dir-ro)-{l-0{^)} 
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V. SUMMARY 



We have shown several numerical results for quantized Duffing oscillator, and discussed 
the crossover from classical to quantum behavior, based on the sensitivity to initial condi- 
tions, " pseudo-Lyapunov exponent ". In our discussion, it is important how A(r) behaves 
as the scaling parameter f3 varies. We have roughly defined (3 = 0.01 as the classical region, 
P = 1.00 as the quantum region and /5 > 1.00 as the deep quantum region. We have found 
the following points, analyzing the numerical results. In the classical region, we prove that 
the chaotic behavior appears certainly, since A(r) increases exponentially. In the quantum 
region, we show that it has been lost completely, since A(r) does not increases exponen- 
tially and takes a certain constant value asymptotically. Notice that there is a clear crossover 
behavior as (3 increases from 0.01 to 1.00; A(r) for (3 = 0.10 still increases exponentially, 
but such a behavior has been lost around (3 = 0.40. For intermediate f3 we can insist that 
the system is classical even in the case of (3 = 0.10. In the deep quantum region, A(r) is 
not a suitable quantity to characterize the behavior of system. It is probably necessary to 
investigate other quantities in order to analyze the behavior in the deep quantum region. 
However, A(r) is an important quantity enough to investigate the crossover behavior be- 
tween the classical and quantum regions, since it clarifies the nontrivial crossover behavior. 
Moreover, we have understood why the chaotic behavior has been lost as /5 1.00 in the 
view of dissipative quantum systems; the effect of dissipation suppresses the occurrence of 
chaotic behavior in the quantum region. The effect of Hamiltonian occur after Tasymp- 
The problem of dissipative quantum chaos has many topics of the foundation for quan- 
^0, Q, [n, 27, 28, 29, Is^. Indeed, there are still many problems for the 



tum theory [ 



quantum-classical correspondence. For example, we have to estimate quantitatively how 
much classical trajectory exerts an effect on the quantum system. A(r) takes a meaning 
value in the classical region, while it is not clear what significant the value of A(r) has in the 
quantum region. Nevertheless, this analysis based on "pseudo-Lyapunov exponent" clarifies 
the crossover in the quantized Duffing oscillator at ffist. This point is quite different from 
the previous workj^. 
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APPENDIX A: THE DERIVATION EQUATION (tTUIl 

We show the derivation of Eq. (jlUj) . To work out this, we should notice that A(r) 
is smaller than the size of effective Planck cell for t > tq. Therefore, it is suggested that 
two wave packets starting from two different initial conditions locate very closely and almost 
overlap respectively. Thus A(r) for r > tq can be approximated by D{Q, P, r) characterized 
by the spread of wave packet for only one initial condition in the phase space. Of course, 
this approximation is valid for r > tq. 



We define {D{Q, P, r)}^ as M{a(a^a)} = M{{a^a) - {a^){a)} with a = [Q + iP) /V2. 



In jl8|, M{a{a\d)} is introduced in order to investigate the localization in framework of 
QSD. Notice that this quantity is different from Trjd^ap} — Tr{a^p}Tr{ap}. M{cr(a^,a)} 
can be only calculated by means of QSD. Considering the average over the set of initial 
conditions, we can write 



fi{Q, P) is the distribution function of initial condition. Since the classical orbits move over 
some finite region in the phase space, e.g., strange attractor, the integration is limited by 
some finite volume V. If fi{Q, P) is uniform distribution, 




(Al) 




(A2) 



If Hamiltonian H = 0, we can derive the following equation: 

-^M{a{a\ a)} = -2TM{a{a\ a)} - 2TM{a{a^ , a)^} 
dr 

-2TM{a{d\d^)a{d,d)} 
< -2^M{(T(a^a)} -2^(M{a(a^a)})2. 
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In the case of equality in the above equations, this relation is Ricaci's differential equation. 
Put u = M{a{a\ a)}, and notice that u is positive for all r. Thus, 



where C = 1 + 1/u{tq). Therefore we obtain the upper bound of A(r) for for r > tq in the 
following equations: 



where we use A(ro) = (3. Eq. (|X3|l is just Eq. (ITn|l . 
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